Sparsity-Aware STAP Algorithms Using Li-norm 
Regularization For Radar Systems 

Zhaocheng Yang, Rodrigo C. de Lamare, Senior Member, IEEE and Xiang Li, Member IEEE 



o 

Oh 






> 

00 

en 

o 

en 



X 



Abstract — This article proposes novel sparsity-aware space- 
time adaptive processing (SA-STAP) algorithms with /i-norm 
regularization for airborne phased-array radar applications. 
The proposed SA-STAP algorithms suppose that a number of 
samples of the full-rank STAP data cube are not meaningful 
for processing and the optimal full-rank STAP filter weight 
vector is sparse, or nearly sparse. The core idea of the proposed 
method is imposing a sparse regularization (/i-norm type) to 
the minimum variance (MV) STAP cost function. Under some 
reasonable assumptions, we firstly propose a /i -based sample 
matrix inversion (SMI) to compute the optimal filter weight vec- 
tor. However, it is impractical due to its matrix inversion, which 
requires a high computational cost when in a large phased-array 
antenna. Then, we devise lower complexity algorithms based on 
conjugate gradient (CG) techniques. A computational complexity 
comparison with the existing algorithms and an analysis of the 
proposed algorithms are conducted. Simulation results with both 
simulated and the Mountain Top data demonstrate that fast 
signal-to-interference-plus-noise-ratio (SINR) convergence and 
good performance of the proposed algorithms are achieved. 

Index Terms — h regularization, Sparsity-aware Space-time 
adaptive processing. Conjugate gradient techniques. Airborne 
radar. Mountain Top data. 



I. Introduction 

Space-time adaptive processing (STAP) is an efficient tool 
for detection of slow targets by airborne or spacebome radar 
systems in serious environments, such as strong clutter and 
lots of jammers [Tl-f4l. However, the full-rank adaptive STAP 
based on linearly constrained minimum variance (LCMV) 
criterion gives rise to two of the major limitations in practical 
applications of radar [T], f4l. First, the computational load 
required to solve the interference matrix inversion is quite 
intense. In addition, the number of training data samples 
required for an accurate estimate of the interference covariance 
matrix can become impractical for high-dimensional problems, 
particularly in heterogeneous environments. It is therefore 
desirable to develop STAP techniques with low computational 
complexity and that can provide high performance in small- 
sample support situations. 

The diagonal loading sample matrix inversion (LSMI) tech- 
nique is considered to be a simple and robust approach for both 
homogeneous and heterogeneous environments ||5], but has a 
high computational cost. Reduced-rank techniques have been 
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investigated for solving the previously discussed problems in 
the last decades |6|-|12|, |14|-|16|, |20|-|22|. One of the 
most important reduced-rank techniques is the class of the 
Krylov subspace methods, which includes the auxiliary-vector 
filters (AVF) |7|, |8|, the multistage Wiener filter (MWF) [9]- 
IT2I and the conjugate gradient (CG) algorithm lfT4l-llT6l. 
1201 . These methods project the observation data onto a lower- 
dimensional Krylov subspace and can obtain an improved 
convergence and tracking performance. The main differences 
amongst them lie on the computational cost, structure of adap- 
tation and ease of implementation. Knowledge-aided (KA) 
STAP techniques have currently gained significant attention as 
an effective STAP algorithm to mitigate the effects of the het- 
erogeneity in the secondary data by exploiting a priori knowl- 
edge II2TI - II24I . However, the exact form of prior knowledge is 
still problem-dependent and hard to be derived. More recently, 
several authors have considered sparse recovery (SR) ideas for 
moving target indication (MTI) and STAP problems Il25l - ||3T| . 
These work based on SR techniques relys on the recovery of 
the clutter power in the angle-Doppler plane, which is usually 
carried out via two steps: first, recovering the clutter angle- 
Doppler profile by some SR algorithms; second, estimating 
the covariance matrix based on the result obtained in the 
first step, and computing the Capon's optimal filter Although 
some fast sparse recovery algorithms are proposed, e.g., the 
fast iterated shrinkage/thresholding (FISTA) algorithm ll27l . 
and the focal underdetermined system solution (FOCUSS) 
based algorithm fSQ\, it is more computationally expensive 
than conventional STAP because of the Capon's optimal filter 
requiring matrix inversion, and the recovery procedure being 
an additive computational burden. 

In airborne radar systems, most interference suppression 
problems are rank deficient in nature |l2]-|0], that is they 
require less adaptive degrees of freedom (DOF) than the full 
DOF provided by the array. In this case, the total adaptive 
DOF provided by the array will be much great than the 
number that needed to suppress the interference. Motivated 
by this, the authors in [32] proposed an sequential approach 
that gave a sparse solution for the transformation matrix to 
select the "best" DOF to be retained in a partially adaptive 
beamformer Moreover, the property described above can be 
seen that there is a high degree of sparsity of the filter weight 
vector Hence, in our prior work, an li type regularization 
to the generalized sidelobe canceler (GSC) STAP processor 
using the ^i-based online coordinate gradient (OCD) method 
[|33]| and the Zi-based recursive least squares method [34] is 
introduced to exploit the sparsity of the received data and filter 
weights, resulting in an improvement in both convergence rate 



and steady-state signal-to-interference-plus-noise ratio (SINR) 
performance. In this paper, we extend the work presented 
in [33J and [jMi to the direct filter STAP processor (DFP). 
By adding the sparsity constraint (Zi-norm regularization) to 
the MV cost function, we derive the Zi -regularized optimal 
filter weight vector under some reasonable assumptions, and 
then propose a sparsity-aware (SA) adaptive STAP strategy 
for airborne radar systems. One direct way is to use the li- 
based SMI recursion algorithm to compute the filter weights. 
However, it requires the matrix inversion operation, which 
prevents its use in practice. The CG method has a low 
computational complexity and is the simplest Krylov subspace 
method since it only needs the forward stage, unlike the MWF 
that requires both forward and backward stages. Therefore, 
low complexity ^i-based CG type algorithms are devised. The 
simulations are conducted using both simulated and measured 
data, which show that the proposed algorithms exhibit im- 
proved performance as compared to existing techniques. 

This paper is organized as follows. Section II introduces 
the STAP signal model for airborne radar In Section III, we 
first introduce the strategy of the SA-STAP algorithm. Then li- 
based SMI and Zi -based CG type algorithms are developed and 
their computational complexity is also shown. Furthermore, 
we conduct an analysis of the proposed algorithms. In Section 
IV, some examples of performance of the proposed algorithms 
with both simulated and the Mountain Top data are exhibited. 
Finally, the conclusions are given in Section V. 

Notation: In this paper, scalar quantities are denoted with 
italic typeface. Lowercase boldface quantities denote vectors 
and uppercase boldface quantities denote matrices. The op- 
erations of transposition, complex conjugation, and conjugate 
transposition are denoted by superscripts T, *, and H, re- 
spectively. The symbols (g) represents the Kronecker product 
and denotes the Hadamard matrix product. Finally, the 
symbol E {■} denotes the expected value of a random quantity, 
operator 3fi[] selects the real part of argument, and the symbol 
Ij • lip denotes the Zp-norm operation of a vector 

II. Signal Model and Problem Statement 

The system under consideration is a pulsed Doppler radar 
residing on an airborne platform. The radar antenna is a 
uniformly linear spaced array (ULA) which consists of M 
elements. The platform is at altitude hp and moving with 
constant velocity Vp. The chosen coordinate system is shown 
in Fig.??. The angle variables (j> and 9 refer to elevation and 
azimuth. The radar transmits a coherent burst of pulses at 
a constant pulse repetition frequency (PRF) fr — 1/T,. , 
where T,. is the pulse repetition interval (PRI). The transmitter 
carrier frequency is fc = c/Xc, where c is the propagation 
velocity and Ac is the wavelength. The coherent processing 
interval (CPI) length is equal to NTr. For each PRI, K 
time samples are collected to cover the range interval. After 
matched filtering to the radar returns from each pulse, the 
received data set for one CPI comprises KNM complex 
baseband samples, which is referred to as the radar datacube 
shown in Fig.??. The data are then processed at one range 
of interest, which corresponds to a slice of the CPI datacube. 



The slice is an M x N matrix which consists of M x 1 spatial 
snapshots for pulses at the range of interest. It is convenient to 
stack the matrix column-wise to form the NM x 1 vector x[A:], 
termed a space-time snapshot, where k is the range sample 
index wdl<k<K fT\-f^. 

Target detection in airborne radar systems can be formulated 
into a binary hypothesis problem, where the hypothesis Hq 
corresponds to target absence and the hypothesis Hi corre- 
sponds to target presence, given as 



Hq : X =Xu 
Hi : X =asS 



(1) 



where as is a complex gain and the vector s, which is the 
NM X 1 normalized space-time steering vector in the space- 
time look-direction, defined as 



s = 



Stifd)®Ss{fs) 
|St(/d)®Ss(/,)|| 



(2) 



where St(/d) denotes the A^ x 1 temporal steering vector 
at the target Doppler frequency fd and Ss{ fs) denotes the 
M X 1 spatial steering vector in the direction provided by the 
target frequency fs- The vector x„ encompasses any undesired 
interference or noise component of the data including clutter 
Xc, jamming Xj and thermal noise x„. Generally, we assume 
the thermal noise is spatially and temporally uncorrelated, and 
the jamming is temporally uncorrelated but spatially strongly 
correlated. As for the clutter, a general model for the clutter 
space-time snapshot is given by l,35J 



Nr Na 



^c[k] = X! X! ^c/fc;m,« (at(fc; m, n) St{ fd/k-7n,n)) 
m—1 71—1 

{as{k; m, n) Ss{fs/k;m,n)) , 



(3) 



where Nr is the number of range ambiguities, Nc is the num- 
ber of independent clutter patches that are evenly distributed 
in azimuth about the radar, (y.t{k; m, n) is a vector describing 
the normalized pulse-to-pulse voltages, and as{k;m,n) ac- 
counts for spatial decorrelation. (7c/k;m.n describes the average 
voltage for the ninth clutter patch and fcth range. The clutter- 
jammer-noise (for short, calling interference in the following 
part) covariance matrix R can be expressed as 



R = ^ {x„x,f } = Re + Rj + R„ 



(4) 



where Re — E{xcX^}, Rj = E{xjXj^] and R„ — 
E |x„x^}, denote clutter, jammer and thermal noise covari- 
ance matrix, respectively. 

Generally, the space-time processor linearly combines the 
elements of the data snapshot, yielding the scalar output [4| 



y = w"x, 



(5) 



where w is the NM x 1 weight vector. The idea behind 
LCMV approach is to minimize the STAP output power whilst 
constraining the gain in the direction of the desired signal. This 
leads to the following power minimization with constraints 



minJ(w) ==^{||w^x||^} s.t. 



w^s 



1. 



(6) 



Using the method of Lagrange multipliers, the optimal full- 
rank LCMV STAP weights are given by ji] 



R 



wlcmv 



ffR-U 



s«R 



(7) 



III. SA-STAP WITH Li-NORM Regularization 

In this section, we detail the design of the proposed SA- 
STAP strategy, derive the /i-based SMI recursion algorithm 
and the Zi-based CG type algorithms and detail their complex- 
ity. Finally, the analysis of the proposed SA-STAP algorithms 
is shown. 

A. SA-STAP Strategy 

In airborne radar systems, most interference suppression 
problems are rank deficient in nature, that is they require less 
adaptive DOF than are offered by the array, the additional 
DOF that are not required can be discarded so that only 
those that are important are retained, which is termed as 
partially STAP technique fSTl. Furthermore, full DOF will 
lead to slow convergence, i.e. requiring many snapshots to 
training the filter, which is difficult to obtain especially in 
non-homogeneous clutter environments. As a result, the total 
adaptive DOF provided by the array will be much great than 
the number that needed to suppress the interference. In another 
word, there is a high degree of sparsity of the filter weight 
vector. However, in the practical, it is not easy to estimate 
the required DOF, related with the sparsity, and to decide 
which DOF are the most important ones. Herein, the authors 
in ll33l . Il34l proposed an li regularized STAP algorithm for 
GSC structure to exploit the sparsity of the received data and 
filter weights. In this paper, we extend this work to a more 
general framework for airborne radar systems, by employing 
the sparse regularization to the MV STAP cost function, which 
is described as the following optimization problem 



mini? • 



{||w^x||2}+2Ar(w) 



s.t. 



w^s = 1, 



(8) 



where A is a positive scalar which provides a trade-off between 
the sparsity and the output interference power The larger 
the chosen A, the more components are shrunk to zero ll36l . 
The sparse regularization is usually conducted by the Zo-norm 
constraint f37l-f39l. However, since this kind of optimization 
problem is known to be NP-hard, one of the approximation 
algorithms, called /i-norm, is considered for the convexity and 
simple complexity 1381 . In the following, we adopt the li- 
norm regularization, i.e., F (v^r) = ||w||i. Now, the question 
that arises is how to effectively solve the li regularized 
MV STAP. Albeit convex, the cost function Ji(w) is non- 
differentiable which leads to difficulty with the use of the 
method of Lagrange multipliers directly. Thus, we propose 
an approximation to the regularization term, which is given 
by 



r(w) 



W 1 



H 



A^ 



where 



A = diag 



|wi| +e' |w2| +e' ' I^TVAfl + e 



(9) 



(10) 



where e is a small positive constant (e.g., e — 0.01 is 
acceptable), and Wi,i = 1,2- •• ,MN are the entries of the 
filter weight vector w. Thus the regularization term w^Aw 
has a quadratic structure, if we assume the diagonal matrix A 
to be fixed. Minimization can be done iteratively by assuming 
that the term A is fixed, being computed with the current 
solution w ll37l . So, fixing the term A, we take the differential 
term with respect to w* of Q, which is given as follows 



5||w| 



,H 



Aw) 



Aw. 



(11) 



9w* (9w* 

The above constrained optimization problem described by 
^ can be transformed into an unconstrained optimization 
problem by the method of Lagrange multipliers, whose cost 
function becomes 

|2" 



C = S {llw^xjl^} + 2A||w||i + 2^, 



fw^s 



1)}, (12) 



where k is a complex Lagrange multiplier Computing the 
gradient terms of (fTZt with respect to w* and k*, we get 



(13) 



V£w* = R-w + AAw + K*s 

V/:„. = w^s - 1. 

By equating the above gradient terms to zero, we obtain the 
filter weight vector 

-1 



w = 



[R + AA] s 



=ff 



R 



AA]"^s' 



(14) 



By inspecting ( fT4b . we verify that there is an additional 
term AA in the inverse of the interference covariance matrix 
R, which is due to the /i-norm regularization. One should note 
that the filter weight vector expression in (fl4l is not a closed- 
form solution since A is a function of w. Thus it is necessary 
to develop an iterative procedure to compute the filter weight 
vector, which will be shown in the following parts. 

B. Li-Based SMI Recursion Algorithm 

In practice, because the interference covariance is unknown 
to us, it is most common to compute the interference covari- 
ance matrix estimate as ll2l-||4l 

1 ^ 

^=iE^Mx''N, (15) 

where {x[n]}^^j^ are known as the secondary or training data. 
In our following derivation, to develop an iterative procedure, 
we add an exponential weighting factor to the interference 
covariance matrix, which may allow the STAP algorithms to 
accommodate possible non-stationarities in the input. We write 
the R[fc] as 



R[fc] 



:/3 



'^-"xfnlx^ln 



[n]x^[n]=/3R[fc-l]+x[i]x^[i], (16) 



where /3 is the forgetting factor, and R[0] = 81, where 5 is 
a small positive quantity and I is the identity matrix. Since 
A.[k] is a function of w[/c], we assume that the filter weight 
values do not change significantly in a single snapshot step, 
which is reasonable because we want the instantaneous error 



of the filter weight vector to change slowly ||40l . Hence, A[fc] 
can be approximated by 

A[k] «A[fc - 1] = 

.. [ 1 1 1 (17) 

However, we note that the computational complexity of 
the ^i-based SMI recursion algorithm is proportional to 
O ( [NM) J , which is not practical, especially in large size 
of phased-array antenna. In the next section, we will develop 
some low complexity algorithms. 

C. Li-Based CG Algorithms 

In order to reduce the computational complexity of the li- 
based SMI recursion algorithm, we introduce low complexity 
adaptive algorithms based on CG techniques to iteratively 
compute the filter weights. There are two different basic strate- 
gies for using the CG method. One is the conventional CG 
(CCG) IfTSl . |fT6l , which executes several iterations per sample 
and runs the reset periodically for convergence. The other 
is the modified CG (MCG) lEl, Ull, ||20l, which operates 
only one iteration per sample. CCG has a faster convergence 
than MCG, but a higher computational complexity. In the 
following, we detail the derivation of the Zi -based SA-STAP 
algorithms based on these two strategies, called ^i-based CCG 
algorithm and ^i-based MCG algorithm. For simplicity, we 
firstly introduce an auxiliary vector given by 



vW 



RL 



AA[fc] 



!■*■ 



(18) 



Then the STAP filter weight vector can be described as w[fc] — 
v[fc]/ (sv[fc]). The solution of v[k] described by ( fTSl l is also 
the solution of the following minimal optimization problem 

mil: 



min JT" (v) — v 



H 



R + AA 



v-2K{v^s}. (19) 



Then the CG-based weight vector is expressed by 

v[fc] = v[k - 1] + a[k]p[k], 



(20) 



where p[fc] is the direction vector, a[k] is the corresponding 
adaptive step size. 

For the ?i -based CCG algorithm, the iteration procedure 
for the CG-based weight vector v is executed per sample. 
For the fcth sample, it assumes constant R[A:] + AA[/c] within 
the internal iterations, and D internal iterations are performed 
per input data sample. The main difference between the li- 
based CCG algorithm and the existing CCG algorithm after 
the derivation is that we add an additional term AA[fc] to the 
estimated interference covariance matrix R[fc]. A summary of 
the algorithm is shown in Table J] 

The /i-based CCG algorithm operates multiple iterations 
per sample and runs the reset periodically for convergence, 
which increases the computational load in the sample-by- 
sample update. In the following, we detail the derivations of 
the Zi-based MCG algorithm with one iteration per sample for 
STAP. From lfT4l . one way to realize the conjugate gradient 
method with one iteration per snapshot is the application 



of the degenerated scheme, which means that the residual 
vector g[k] will not be completely orthogonal to the subspace 
spanned by the direction vectors {p[0],p[l], • ■ • ,p[fc — 1]}. 
Under this condition, the adaptive step size a[k] has to fulfill 
the convergence bound given by 

< |p^[fc]g[A:]| < 0.5 |p^[fc]g[fc - 1]| , (21) 

where g[fc] is the negative gradient vector of J' (v) in iT% . 
Thus, g[fc] can be written as 



g[fc] = -V J (v)^, = - [R[fc] + AA[fc] 
which can be calculated recursively by 

g[fc] = (1 - /3) s + /3g[fc - 1] - a[k] \R.[k] 



r[k] 



(22) 



AA[fc] p[fc] 



(23) 



- {A [1 - 13] A[k] + x[fc]x^[A:]} v[fc - 1]. 

In the previous equation, we use the approximation that 
A[fc — 1] w A[fc]. Premultiplying (l23T l by p^[fc], taking the 
expectation of both sides and considering p[fc] uncorrected 
with s, x[fc] and v[/c — 1] |14J, we obtain 

E[p"[k]g[k]] ^l3E[p"[k]g[k-l]] 

-/3E[p"[k]s] (24) 

- E \a[k]p"[k] fR[fc] + AA[A:]) p[k] 



Here, it is assumed that the algorithm converges 
with the assumption that E [v[fc — 1] — Vopt] ~ 0, 
E [x[fc]r^[/c]v[fc - 1]] « s, and E[X[1- /?] A[k]^r[k - 1]] « 
0. Making a rearrangement of (l25t and following the 
convergence bound (l2Tl i. we obtain 

'1 



a[k] 



p"[k](R[k]+XA[k]\p[k] 
-p"[k]s]-^ip"[kMk-l]}. 



{P[p"[k]g[k-l] 



(25) 



where < fi < 0.5. The direction vector is a linear com- 
bination from the previous direction vector and the negative 
gradient, which is described as 



p[k]=g[k~l] + iy[k]p[k], 



(26) 



where i^[fc] is computed for avoiding the reset procedure by 
employing the Polak-Ribiere approach, which should have an 
improved performance llT4l . llT6l . and is stated as 



H 



v[k] 



[g[k] - g[k - 1]]" g[k] 



(27) 



gfi[k-l]g[k-l] • 

The proposed /i-based MCG STAP algorithm is summarized 
in Table HIl 

From above discussions, two aspects should be noted that: 
First, the performance of our proposed algorithms (both li- 
based SMI and /i-based CG-type algorithms) depends on 
regularization parameter A. An approach to choose A is 
introduced in ll34l . which can be easy to extend to our 
proposed algorithms, but not discussed in this paper for saving 
space. Second, the convergence analysis in |16| is suitable to 
our proposed CG-type algorithms, where the convergence is 
governed by 



K»+i| 



lG[fe] 



<2 



V '^max/'^i 



max/ ' inin 



/Tit 



ll^oi 



|G[fc], (28) 



' max/ ' mm 



where ^^ [k] = Vopt [k] — v^ [k] is the CG-based weight vector 
error at the ith iteration for the fcth snapshot, Vopt[fc] is the 
optimal solution at the fcth snapshot, Tj^hx and Ty^in are the 
maximal and minimal eigenvalues with respect to G[fc] = 
R[fc] + AA[fc], and ||Q[fc]|lG[fc] = <;f[k]G[k](;,[k]. From the 
above equation, we note that the convergence behavior of the 
proposed algorithms is related to the CG-based weight vector 
error ^o[^] and the condition number Tj^hx/tu 



'max/ 'nun* 



D. Complexity Analysis 

In this section, we detail the computational complexity in 
terms of complex additions and complex multiplications of 
the proposed ^i-based SMI, /i-based CG type algorithms, and 
other existing STAP algorithms, namely the LSMI, the AVF, 
the MWF and the conventional CG type algorithms, as shown 
in Table |III1 One aspect should be noted that, the rank D may 
not equal to the clutter rank, and can be smaller than that. 
This is because the principle of the Krylov subspace approach 
is different from that of the eigen-decomposition approach. 
An eigen-decomposition approach would usually require an 
SVD on the full-rank covariance matrix and the selection of 
the D eigenvectors associated with the D largest eigenvalues, 
which is high related to the clutter rank. In contrast to that, the 
Krylov-based approach does not require eigen-decomposition 
and selects the D basis vectors which minimize the desired 
cost function and will form the projection matrix, where D 
can be decreased without significantly degrading the SINR 
ll^. In the table, D is the rank for CCG type, AVF and 
MWF algorithms, and L = NAI is the system size. Seen 
from the table, the computational complexity of Zi -based SMI 
is similar to the conventional LSMI algorithm, both requiring 
one to calculate the matrix inversion. With respect to the 
proposed Zi-based CG type algorithms, the computational 
complexity is nearly the same as the conventional CG type 
algorithms. Note that the complexity of CCG type, AVF and 
MWF algorithms is dependent on the rank D. This is a 
tradeoff between complexity and performance. We found that 
the rank of the proposed ^i-based CCG algorithm with D = 7 
works well (while the best rank for AVF and MWF is much 
larger), as will be verified in the following simulations. The 
low-rank characteristic will bring computational savings. The 
computational complexity of all algorithms is shown in FigUl 
where we use the best rank obtained from the simulations for 
these algorithms (D = 7 for CCG type, D = 18 for AVF 
and D = 14 for MWF). We see that the proposed CG type 
algorithms have much lower complexity than AVF and MWF 
algorithms. 

Furthermore, it requires to compute the filter weights repeat- 
edly for target detection in airborne radar systems, especially 
in heterogeneous environment. In this case, our proposed 
algorithms can work in an iterative way and do not need to 
recompute all the filter weights, which can lead to significant 
computational savings. Usually, secondary data of the sliding 
window are used in detection procedures, where the parameter 
that defines the length of the sliding window is K. Assume 
Ri[iir] denotes the estimated interference covariance matrix 
according to (fTSI l and Wi[A'] denotes the filter weight vector 



at the cell under test (CUT) of the ith range bin, respectively. 
Consider the case of the i + 1th CUT, we first remove the 
impact of i + 1th CUT, given by 



/?R,+i[if-l]=R,[A']-x[i- 



l]x«[ 



(29) 



Since an exponentially decaying data window is used, we do 
not need to remove the first snapshot used to compute the 
filter weights. Then, similarly, we consider the case of adding 
snapshots. Two snapshots, one is at the primary ith CUT and 
another is the new snapshot x„cw which was not included in 
the sliding window before, should be added to the i + 1th CUT 
secondary data. The procedure can be written as 



R,+i[K] ^(5Ii,+i[K - l]- 



/3x„owxfew + x[i + l]x"[i + 1]. 



(30) 



As for the filter weight vector Wj+i [K] at the i + 1th CUT, it 
can be updated using the new interference covariance matrix 
and the filter weight vector w^ [X — 1] . 

In addition, the proposed algorithms adopt an adaptive filter- 
ing approaches, which can obtain a near optimum interference 
rejection at a low cost [4iJ. The advantage of this approach 
is that filtering can be accomplished in a pipeline mode as 
the echo pulses come in. The required number of calculations 
for filtering can be realized easily with nowadays digital 
technology ll42l . 



E. Analysis of the SA-STAP Algorithm 

At this point, we have finished the derivation of the SA- 
STAP algorithms. The following simulation results will show 
that the proposed SA-STAP algorithms have a faster SINR 
convergence speed and better SINR steady-state performance 
than the conventional algorithms. This translates into a supe- 
rior detection performance. However, why do the SA-STAP 
algorithms work is an interesting question. This section will 
try to explain that from two points of view. 

First, to understand the behavior, we write the filter weight 
vector using the eigenvalue decomposition (EVD) of R. We 
assume that the eigenvalues of the estimated interference 
covariance matrix are 7„ with the corresponding eigenvectors 
denoted by u„, n = 1,2,-- ,NM. The eigenvalues are 
ordered as. 



71 > 72 > ■ ■ • > INM = In 



(31) 



Thus, through the EVD, the estimated interference covariance 
matrix can be written as 



NM 

R= ^7„u„uf. 



(32) 



Substituting ( [32] i into ( fT4b . the filter weight vector of the SA 
STAP algorithm can be written as 

( NM 



'u„ rt) u„ \ , (33) 



fl % + Sit + ^n 

Where S^t = miii(^^),ri = 1,2,---,NM, A. 

the difference between -, — 4-;— and (5£^„, and qsA is a sc 
quantity, which does not affect the SINR. 



IS 



By inspecting ( [33T l. we observe that the SA-STAP belongs 
to the class of diagonal loading STAP techniques in a sense. 
Moreover, it is equivalent to an adaptive diagonal loading tech- 
nique, which will apply to each eigenbeam of the interference 
CO variance matrix different weights and exploit the sparsity of 
the filter weights and the received data. 

Second, we will investigate the relationship between the 
SINR performance and the Zi -norm-sum quantity of the filter 
weights. Assume the scene is the same as the one with 
homogeneous environment introduced in the next section. We 
compute the SINR loss and the ^i-norm-sum quantity of the 
filter weights against the number of snapshots using the SMI 
algorithm. The results are plotted in FigjS] From the figure, 
we find that the better the SINR performance, the smaller the 
Zi-norm-sum quantity of the filter weights. From this point of 
view, a constraint on the /i-norm-sum quantity of the filter 
weights can lead to fast convergence, which in fact exploits 
the sparsity of the received data and filter weights. 

IV. Performance Assessment 

In this section, we assess the proposed SA-STAP algorithms 
using both simulated and measured data and compare them 
with the existing algorithms, such as the conventional CCG, 
MCG, MWF AVF and LSMI algorithms. We measure the 
SINR, the SINR loss and the probability of detection curves, 
where the SINR and the SINR loss are defined as follows |4], 
respectively. 

I -ff |2 
W X 

SINR = ,L ,._L, . (34) 



tH 



Rwl 



rH^ 



SINRlosS = 



rH 



Rwl 



(35) 



where R is the exact interference covariance matrix at the 
detection range bin and w is the estimated filter weights using 
the neighbor secondary data. 

A. Simulated Data 

Consider a monostatic sidelooking radar with M = 10 
antenna elements and A^ = 8 pulses in one CPI, giving a 
space-time steering vector of length L — 80. We assume 
a simulated scenario with the following parameters: half- 
wavelength spaced antennas, uniform transmit pattern, carrier 
frequency 45GMHz, PRF set to 300Hz, platform velocity of 
50m/s and height of 9000m, the clutter uniformly distributed 
from azimuth — 7r/2 to 7r/2 with clutter-to-noise-ratio (CNR) 
of 40dB, two jammer located at —45 and 60 with jammer-to- 
noise-ratio (JNR) of 40dB, the target located at 0° azimuth 
with Doppler frequency of lOOHz and signal-to-noise-ratio 
(SNR) of OdB, and the thermal noise power is O.OIW. We 
consider the inner clutter motion (ICM) in simulated data. 
One common model, referred to as the Billingsley model, was 
developed by Billingsley of MIT Lincoln Laboratory |3|. The 
only parameters required to specify the clutter Doppler power 
spectrum are essentially the shape parameter b and the wind 
speed parameter w. In this paper, we assume h — 3.8 and 



Lo = 51.45 miles per hour (mph). All presented results are 
averages over 100 independent Monte Carlo runs. 

In our first example, we consider the SINR performance 
versus the rank D of the proposed /i-based CCG algorithm, 
the conventional CCG algorithm, the AVF algorithm and the 
MWF algorithm. A total of K — 160 snapshots are considered. 
The results in Fig|3] show that our proposed Zi-based CCG 
algorithm can obtain its best performance when the rank is 
larger than _D = 7. It is much lower rank to obtain its best 
performance than that of AVF {D = 18) and MWF {D = 14) 
algorithms. The low-rank characteristic will bring considerable 
computational savings, which is very important for STAP in 
radar systems. One should note that the performance of the 
conventional CCG algorithm will degrade when the rank is too 
large, while our proposed ?i -based CCG can always keep good 
performance resulting in further robustness. Since the SINR 
performance is much worse when the rank is lower than the 
best rank, thus, we will use D = 7 for CCG type algorithms, 
D = 18 for the AVF algorithm and D = U for the MWF 
algorithm in the following examples. 

In the next example, we evaluate the SINR loss performance 
against the number of snapshots K = 320 of the proposed 
algorithms with the existing algorithms, as depicted in Fig|4] 
The curves show that: (1) the SINR performance of the pro- 
posed /i-based SMI algorithm is a suboptimal algorithm, but 
exhibits the best performance compared with other algorithms. 
(2) Zi-based CG type algorithms outperform conventional CG 
type algorithms in terms of convergence rate and steady- 
state performance; (3) the SINR performance of the /i-based 
CCG algorithm is better than AVF and MWF algorithms. 
(4) Although the Zi -based MCG algorithm shows slower 
SINR convergence than the MWF algorithm, we can obtain 
a better SINR performance when the number of snapshots is 
larger than 100. One should note that the proposed CG type 
algorithms have a much lower computational complexity than 
LSMI, AVF and MWF algorithms. 

In the third example, we present the probability of detection 
Pd versus SNR with the target injected at the azimuth of 0° 
and Doppler frequency lOOHz in Fig|5] We assume the false 
alarm rate Pfa is set to 10^^ and the number of the secondary 
data is K — 110. The plots illustrate a similar trend to the 
SINR loss performance in the second example. Note that we 
obtain a performance gain of about IdB in terms of SNR for 
Zi-based CG type algorithms, as compared with conventional 
CG type algorithms. 

Fig|6] shows the SINR performance against the target 
Doppler frequency at the azimuth of 0° with a total of 
K = 160 snapshots. Here, we suppose the potential Doppler 
frequency space is from — lOOHz to lOOHz. The parameters 
of all algorithms are the same as the second example. The 
curves in the figure demonstrate a similar trend to the results of 
previous examples. Additionally, the Zi-based SMI algorithm 
displays much better performance to the slow targets than other 
algorithms. 

B. Measured Data 

In this section, we apply the proposed algorithms to the 
Mountain-Top data set. This data was collected from com- 



manding sites (mountain tops) and radar motion is emulated 
using a technique developed at Lincoln Laboratories f6], ['431. 
The sensor consists of 14 elements and the data are organized 
in CPIs of 16 pulses. Here, we use the data file t38pre01vl 
CPI6, which could be obtained from the internet (A4]. The 
pulse PRF was 625Hz and the instance bandwidth after pulse 
compression was 500kHz. There are 403 independent range 
samples available for the training data support. The clutter 
was located around 245° azimuth and the target was at 275° 
with a Doppler frequency 156Hz. All the data processed 
following are through pulse compression firstly. Note that the 
clutter and target have the same Doppler frequency, hence 
separation is impossible in the Doppler domain but possible in 
the spatial domain. The estimated angle-Doppler profile using 
all 403 samples is given in Fig.??, which shows a serious 
heterogeneity. 

Fig. and Fig. display the STAP output power of all algo- 
rithms in the range of 147-162 km. Here, the interference 
covariance matrix estimated using a symmetric sliding window 
with a total of 20 snapshots for Fig. and 40 snapshots for Fig.. 
For each CUT, the snapshots do not include the 6 snapshots 
around the CUT. In the figures, we also give the unadapted 
weight vector, which equals the steering vector w = s. We 
see that the target is clearly not detectable without adaptive 
processing. To have a clear comparison amongst different 
algorithms, we show the differences between the output power 
at 154 km and the next highest power peak in Table |IV] 
where "-" presents the target not detectable. Here, 6 range bins 
around the rang bin of the target is not used for comparison 
since they are the guide cells. Seen from the table, we find 
that: (1) the proposed ?i -based SMI algorithm obtains the 
best detection performance in both situations, which is the 
same conclusion as that using simulated data; (2) the proposed 
Zi -based CG type algorithms obtain better performance than 
the conventional CG type algorithms (although the proposed 
Zi -based MCG algorithm has a pseudo target at the range 
153km, when the secondary data record is 20 snapshots, the 
conventional MCG algorithm can not detect the target at 
all.); (3) the proposed Zi-based CCG algorithm outperforms 
AVF and MWF algorithms in both situations. Hence, we 
can conclude that our proposed algorithms show a robust 
performance in heterogeneous environments. 

V. Conclusions 

In this paper, we have proposed novel SA-STAP algorithms 
with Zi-norm regularization for targets detection in airborne 
radar systems. The proposed SA-STAP algorithms employed 
a sparse regularization to the MV cost function to exploit the 
sparsity of the received data and filter weights. To solve this 
kind of optimization problem, an /i-based SMI algorithm was 
directly developed, but it required matrix inversion resulting 
in a high computational cost. Accordingly, we have proposed 
low-complexity SA-STAP algorithms based on CG techniques. 
A detailed analysis of the computational complexity and 
the performance of the SA-STAP algorithms were carried 
out. Simulation results with both simulated and measured 
data showed that the proposed algorithms outperformed con- 



ventional algorithms and exhibited a robust performance in 
heterogeneous environments. 
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TABLE I 
The h-BASED CCG Algorithm 



Initialization: 

R[0] = 51, v[0] 



■ s,rjcca. 



Recursion: For eacli snapsliot k = 1, ■ • • , L 
STEP 1: Start: 

R[fc] = /3R[A: - 1] + x[fc]x^[fc], 



G[k] = R.[k] + XA[k], 

go[fc] = s - G[fc]vo[fc], pi[fc] = go[fc], Po[k] = g(f [fc]go[fc], 
STEP 2: For d = 1, • • • , D and Pd-i[k] > r?ccG 

z4k] = G[k]p4k], 

a4k]^ [p^[k]z4k]y'pk-i[k], 
Vd[fc] = Vd-i[fc - 1] + ad[k]pd[k], 
Sdlk] = gd-i[fc] - Qd[fc]zd[fc], 
p4k]^g^[k]g4k], 

-4fc] = tA-v 

Pd+i[k] = gdW + i^d[k]pd[k], 
STEP 3: After end STEP 2 

Vo[fe+l] =Vd[fc], 

w[fc] 



vqW 
s-"vo[fcl' 



Final output: 

y[k] = w^[fc]x[fc]. 



TABLE II 
The /i-based MCG Algorithm 



Initialization: 

RJO] = SI, v[0] 



s, w[0] = s, g[0] = s, p[l] = s, 



,L 



1 \ 

fe-lJI + Ej' 



Recursion: For each snapshot k 

R[fc] = /?R[fc - 1] + x[fc]x^[fc], 

G[k] =R[fc] + AA[fc], 
a[k]=[p"[k]G\t]p[k]]-' 

X {13 [p"[kUk - 1] - p"[k]s] - pp"[k]g[k - 1]}, 
v[fc] = v[j — 1] + a[i]p[i], 
g[fcl = {1-P)s + Ps[k - 1] - a[k]G[k]p[k] 

-~{\[1~ l3]A[k] +x[fc]x"[fc]}v[fc-l], 
,,rw _ [g[fc]-g[fc-i]]^g[fc] 

^ > S"lk-l]g[k-l] ' 

p[k + l] = s[k - I] + 4k]p[k], 
w[k] = -^^^ 



7TT 



vW 



Output: 

y[k] = w^[fc]x[fc]. 



TABLE III 

Comparison of the Computational Complexity 



Algorithm 


Additions 


Multiplications 


LSMI 


O (L^) + O [L') 


O (L^) + O (L') 


h -based SMI 


0(L-^)+0(L^) 


0(L^)+0(L^) 


MWF 


DL'^ + {AD - D'^)L 

,D^ 313^ D 
'' -.^ '>. ?. 


2DL'^ + {5D - D'^)L 
+ 2f 2^2 +16° 


AVF 


(2D + 1)1-' + {AD 
+1)L-AD- 1 


2{D + 1)L-' 
+7DL + L 


CCG 


{D + 2)L-' + {AD 
+2)L -2D -2 


{D + 3)L' 
+5DL + 3L 


MCG 


3L^ + lOL - 4 


AL^ + 13L + 2 


ii -based CCG 


(D + 3)L" + {AD+ 
3)L - 2D - 2 


{D + 3)L' 
+7DL + L 


ii -based MCG 


5L" + llL-4 


AL^ + 13L + 2 



TABLE IV 

Results of Mountain Top data 



Algorithms 


20 snapshots 


40 snapshots 


unadapted 


- 


- 


LSMI 


7.1dB 


12.4dB 


ii -based SMI 


9.3dB 


IS.ldB 


MWF 


4.6dB 


lO.OdB 


AVF 


5.6dB 


11.9dB 


CCG 


5.2dB 


11.9dB 


MCG 


- 


6.2dB 


/i -based CCG 


7.1dB 


IS.OdB 


Zi -based MCG 


OdB 


9.5dB 





Fig. 3. The SINR perfoiTnance versus the rank D. Parameters: the diagonal 
loading factor for AVF and MWF algorithms is lOdB to the thermal noise 
power; /3 = 0.9998, J7CCG = 10"^ and R[0] = O.OOII for CCG type 
algorithms; A = 2 for /i -based CCG algorithms. 
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Fig. 1. The computational complexity per snapshot. 
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Fig. 4. The SINR loss performance against the number of snapshots 

K = 320. Parameters: the diagonal loading factor for LSMI, AVF and 
MWF algorithms is lOdB to the thermal noise power; /3 = 0.9998 and 
R[0] = O.OOII for CG type algorithms; A = 1 for the /i-based SMI, A = 2 
for the /i-based CCG algorithm and A = 1 for the h-based MCG algorithm; 
^CCG = 10"''. 



Fig. 2. The relationship between the SINR performance and the Ji -norm- 
sum quantity of the filter weights. 
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Fig. 5. Probability of detection performance versus SNR with K = 110 
snapshots. Pf^ = 10^'' and the other parameters are the same as the second 
example. 
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Fig. 6. SINR performance against Doppler frequency with K = 160 
snapshots and Doppler frequency space from —100 to lOOHz. The other 
parameters are the same as the second example. 



